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Oxide heterostructures are of great interest both for fundamental and applicative reasons. In 
particular the two-dimensional electron gas at the LaAlOa/SrTiOa or LaTiOa/SrTiOa interfaces 
displays many different physical properties and functionalities. However there are clear indications 
that the interface electronic state is strongly inhomogeneous and therefore it is crucially relevant to 
investigate possible intrinsic electronic mechanisms underlying this inhomogeneity. Here the electro¬ 
static potential confining the electron gas at the interface is calculated self-consistently, finding that 
the electron confinement at the interface may induce phase separation, to avoid a thermodynam¬ 
ically unstable state with a negative compressibility. This provides a generic robust and intrinsic 
mechanism for the experimentally observed inhomogeneous character of these interfaces. 

PACS numbers: 73.20.-r,73.43.Nq,73.21.Fg, 74.81.-g 


The two-dimensional electron gas (2DEG) that 
forms at the interface of two insulating oxides, like 
LaAlOa/SrTiOa and LaTiOa/SrTiOs (hereafter gener- 
ically referred to as LXO/STO) [IHl], exhibits a 
rich phenomenology, such as a gate-tunable metal-to- 
superconductor transition [SHE], a magnetic-field-tuned 
quantum criticality [5] , and inhomogeneous magnetic re¬ 
sponses [lOHIb] . Tunneling [m E] and SQUID mag- 
netometry [18] provide clear evidence of an inhomoge¬ 
neous interface on both micro- and nanoscopic scales. 
Transport measurements report further signs of inhomo¬ 
geneity and a percolative metal-to-superconductor tran¬ 
sition with a sizable fraction of the 2DEG never becoming 
superconducting down to the lowest accessible tempera¬ 
tures [TUIE2] . For both fundamental reasons and applica¬ 
tive purposes, like device design, it is crucial to identify 
possible intrinsic mechanisms that may render the 2DEG 
so strongly inhomogeneous via a phase separation (PS). 
This is precisely the focus of the present work. 

Here, we identify a very effective electron-driven mech¬ 
anism leading to PS, based on the confinement of the 
2DEG at the interface. From customary self-consistent 
calculations of the conhning potential well in semicon¬ 
ductors, it is well known [23] that a hnite lateral exten¬ 
sion usually renders the 2DEG more compressible than 
its strictly 2D counterpart. This effect is much stronger 
in LXO/STO than in ordinary semiconductor interfaces, 
due to the huge dielectric constant of STO, allowing for 
much larger electron densities, with a strong amplifica¬ 
tion of the self-consistent adjustments of the conhning 
potential. As a consequence, a non-rigid band structure 
arises, which varies with the local electron density: an in¬ 
creased electron density is accompanied by a correspond¬ 
ing increase of the positive countercharges (due to oxy¬ 
gen vacancies and/or polarity catastrophe [2^29] ). from 
which the interfacial electrons are introduced and restor¬ 
ing the overall charge neutrality. For small-to-moderate 


increases of electron and countercharge densities the po¬ 
tential well deepens and the electron energy levels are 
shifted downwards. In this Letter we show that this 
mechanism leads by itself to PS. 



FIG. 1. Sketch of the interface for back (a) and top (b) gating. 
The confining potential depends on both mobile (nlf) and 
trapped (uq) charges. Applying a positive (negative) voltage, 
Sn electrons per unit cell are added to (subtracted from) the 
interface and the potential changes accordingly. 

— The model — The thermodynamic stability of the 
system is investigated by varying the density of the inter¬ 
facial gas while keeping the overall neutrality. Therefore 
a corresponding amount of positive countercharges has 
to be varied (see Fig.[^. Due to this tight connection be¬ 
tween positive and negative charges the (in)stability will 
be determined by calculating the chemical potential of 
the whole system (i.e., of both the mobile electrons and 
of the rest of the charges). While for the electron part we 
will solve the quantum problem of the mobile electrons 
in the self-consistent confining well, the countercharges, 
the fraction of electrons trapped in impurity states of the 
bulk (see below) and the boundary conditions fixing the 
gating potential will determine the classical electrostatic 
energy of the system. We will keep all these contribu¬ 
tions (for a detailed description see the Appendix) into 
account to calculate the total energy £ and, in turn, the 
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chemical potential ^ = £(N +1) — £{N) « dM£ (here N 
represents the number of electrons, which is always kept 
equal to the number of countercharges). 

The energy of the 2DEG is obtained through the calcu¬ 
lation of the energy spectra as a function of the electron 
density, n = Sn + no [henceforth, carrier densities, e.g., n, 
carrying no z dependence are meant per interfacial unit 
cell (u.c.), and are related to their z-dependent counter¬ 
parts by relations like, e.g., n = dz n(z)]. Here, the 
two contributions have different origin: Sn is the extrin¬ 
sic component introduced by gating, while the intrinsic 
component uq 1301 originates from the electronic recon¬ 
struction due to the polarity catastrophe and/or from 
oxygen vacancies within the LXO layer. Which of the 
two dominates is not important in our calculations and 
we represent the related positive countercharges as the 
light-blue shaded areas in the LXO side of Fig.[^ 

What matters, instead, is the way the extrinsic charges 
are introduced, particularly in the case of back-gating 
[Fig.[^(a)]: Applying a positive voltage Vg, the electro¬ 
static potential, after increasing in the region close to 
the interface, decreases linearly with distance, once the 
interfacial charge density has been exhausted [the elec¬ 
trons reside on the STO side, which we assume to occupy 
the z > 0 half space and n{z) —)■ 0 for large z]. Then the 
electrons closest to the top of the well are weakly con¬ 
fined and some of them may escape and get trapped into 
the defects of bulk STO, as it is demonstrated by irre¬ 
versibility effects under large back-gating potentials [3T]. 
Thus in the absence of trapped charges, the quantum 
well is intrinsically unstable upon positive back-gating. 
In the top-gating configuration [Fig.[^(b)], the leakage 
also occurs, because the Fermi energy of the electrons, 
attracted to the interface by the positive Vg, exceeds the 
confining potential on the STO side. In both configura¬ 
tions, we are therefore led to introduce trapped charges 
that we describe by a distribution nQ(z) = (ng/A)e“^/'^ 
on the STO side, decaying over a distance A of several 
tens of nanometers (see pink shaded area in Fig.[^. This 
has the main effect of deepening the confining potential 
well, the Fermi energy being located substantially below 
its top (see, e.g.. Refs. 151 and 1^ . 

The mobile electrons occupy energy levels that are 
quantized in the z direction and form a 2D band struc¬ 
ture in the xy (interfacial) plane. The electrostatic po¬ 
tential (j){z) confining the electrons is determined self- 
consistently with the mobile electron density distribu¬ 
tion n‘^{z) = nJJ^(z) + Sn(z) (at external gating Vg), for 
a frozen distribution of trapped charges 'Uq(z). The z 
component of the factorized wave function ^(x,y,z) = 
C(z)'ipk^ky(x,y) is the solution of the Schrodinger equa¬ 
tion yielding the sub-band energy e^. 


■ 

_2m^ dz"^ 


-I- e(t>{z) + Eia 


Cia (^) — 0 ; 


( 1 ) 


where the electron charge is —e and i = 1, 2, 3,... is the 


sub-band index. The index a = xy, xz, yz labels the 
three Ti t 2 g orbitals, d^y, d^z-, dyz, where the electron 
mostly reside. The full 2D band structure is 

_ h^kl , ^ _ 

^2,Q,k — ... T V zAq, -t- Sia- 

2mg znia 

Taking = 0, the energy offset A^z^yz = A of the 
dxz,{yz) bands is experimentally found to be A « SOmeV 
[34) . Similar values are found in first-principle calcula¬ 
tions between the highest occupied dxy sub-band and 
the lowest dxz,yz sub-bands [35]. Here we have to ad¬ 
just A « 0 — 10 meV in order to self-consistently obtain 
the energy difference between the dxy and the d^z^yz sub¬ 
bands of order 50 — 60meV (see the Appendix), showing 
that the energy offset mostly arises from the different 
mass along the z direction. We take the masses of the 
various bands as = to/, m%y = mh', 1 x 1 %’^ = mi, 
m^z = ixih', = TO/, TOg^ = TO/j. According to stan¬ 

dard values [35], we take to/ = 0.7 me and mh = IdTOe 
{me is the electron mass). Assuming full translational in¬ 
variance along the xy planes and integrating over a u.c. 
of area with a suitable normalization of 'ipk„:ky{x,y), 
the density of mobile electrons, at temperature T = 0, 
reads 


'(^) = / deg^aie), 

J —OQ 


where gia (e) = a?{nfi^)~^^/m%ma 0{e — Sia — A^) is the 
density of states (DOS) of the various Og sub-bands, 0{e) 
is the Heaviside function, and Ef is the Fermi energy. 

The electron distribution corresponds to an electro¬ 
static potential ^(z) obeying the Poisson equation: 

= n"^{z) + nl{z). (2) 

Here, the dielectric constant is a function of the elec¬ 
tric field E = —d(j)/dz via the relation er{E) = 
{A + B\E\) -I- eooj where A, B, and Coo are experimen¬ 

tally measured constants m- Owing to the nearly ferro¬ 
electric character of STO, Er can reach very large values 
(> 25x 10^) but, due to the very strong interfacial electric 
field [Sl[35], we find that, near z = 0, ej. « Coo ~ 100—300. 
The crucial point of the above derivation is that the cal¬ 
culation is self-consistent only if the two potentials, (j) 
from Eq. Q and (j) from Eq. ([^, coincide. 

Besides the difficulties stemming from the self- 
consistency, there are additional subtleties coming from 
the boundary conditions, which vary for the back- or top¬ 
gating configuration. In the former case we fix a density 
of positive charges uq = dz [n™(z) -|- Uo(z)] (per u.c.) 
at z = 0“. We made this simplifying choice, which turns 
out to be slightly less favorable to the occurrence of PS, 
to avoid the distinction between the case of oxygen va¬ 
cancies (for which the countercharge would be uniformly 


eoa"' 


dz 


er{E)—(j){z) 

dz 
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distributed in the LXO) and the polarity catastrophe (in 
which the countercharges suitably distribute themselves 
in the polar planes of LXO in order to minimize the en¬ 
ergy [53] )• Having fixed the positive charges, we consider 
the electric field at 2 = O'*' (i.e., the slope of the con¬ 
fining potential) at the interface. The electric field deep 
inside the STO [where no{z) 0] is fully determined by 
the gate potential because the intrinsic electronic density 
and the corresponding positive charges in the LXO com¬ 
pensate and have thus no effect in this region. On the 
other hand, the electrons coming from back-gating create 
a field Vg/L, with eSn = er{E)dE, where L is 

the thickness of the STO substrate. In the top-gating 
configuration, instead, the density of positive charges at 
z = 0~ amounts to no + Sn, and the electric field vanishes 
deep inside the STO substrate. 

Once the quantum problem of mobile electrons is 
solved, the rest of the energy is due to the electrostatic 
contribution (per u.c.) of electric fields due to the over¬ 
all distribution of the mobile (m) and the fixed (/, from 
gates and trapped) charges 

, „2 rca 

= er{E) [Ef^iz) - E^^iz)] dz, (3) 

with E = E^ + E"^. Notice that, since the Hartree- 
like electrostatic energy is double-counted in the quan¬ 
tum Hamiltonian, the contribution of the mobile charges 
must be subtracted in Eq. (|^. In the Appendix we pro¬ 
vide details of how the above fields are calculated. 

— Results — To evidentiate that the electron confine¬ 
ment is the driving mechanism of PS we first consider 
the pure quantum problem neglecting the classical elec¬ 
trostatic contributions to the free energy. We thus report 
in Fig 2(a) the evolution of the different sub-bands levels 
as a function of the mobile electron density n™, which 
is the sum of the initial as grown carrier density n™ and 
the extra charges 6n added by electrostatic gating. This 
evolution is a direct consequence of the non-rigidity of 
the bands. The resulting Fermi energy Ep is going up as 
expected, while the chemical potential p,™ of the mobile 
electrons has a non-monotonic behavior. Indeed, when 
Ep crosses the first heavy band {dxz,yz which have the 
largest DOS along the x — y plane), /r’" starts to de¬ 
crease. The most important result of this paper is that, 
for a given range of mobile electron densities, the chemi¬ 
cal potential decreases while electrons are added, result¬ 
ing into a negative compressibility. The entailed ther¬ 
modynamic instability is prevented by PS, with the con¬ 
sequent inhomogeneous redistribution of electrons and 
countercharges. Remarkably, the electron system tends 
to stay unstable (i.e. it has a negative compressibility) 
up to very large densities showing that this instability 
mechanism is strong and robust. However, this cannot 
be a physical situation, since the compressibility remains 
negative, preventing a stable phase separation to occur. 



FIG. 2. (a) Fermi energy (magenta solid line), electron-only 
(see text) chemical potential (blue solid line), sub-band lev¬ 
els (green dashed for the d^y and black dashed for the dxz,yz 
levels) as a function of the build-in density no at zero gating 
potential, (b) Chemical potential as a function of the mobile 
electron density n™ = -|- 5n at fixed values of the back 

gating potential Vg oc Sn (the Sn electrons due to gating are 
thus also fixed). The black solid (red dashed) lines correspond 
to the solution in the absence (presence) of a countercharge 
background with short-range rigidity. The short-range back¬ 
ground contribution to the chemical potential is taken of the 
form fisr{n^) = A [(n’" — Sn) /0.065]^, where the exponent 
p = 19 measures the short-range rigidity and A = 7 x 10“^. 
For Vg = 50 V, Sn = 0.013668. (c) Same as (b) with Sn = 0. 
(d) Same as (b) with Sn = —0.01711. 


Indeed, this model is not complete. One has to introduce 
the energetic cost of the electrostatic fields generated by 
all the charges in the system according to Eq. (§. This 
leads to the chemical potential displayed in the panels of 
Fig. i (b-d) (black solid line) when n™ is increased (the 
gate voltage is kept constant in this calculation, while the 
intrinsic n™ is varied).. Taking into account the quan¬ 
tum and electrostatic energies restores a physical situa¬ 
tion, where the negative compressibility occurs in a lim¬ 
ited range of doping, since for high n™, y increases when 
adding carriers. 

Now the increasing classic electrostatic energy cost 
tend to stabilize the system at large and allow to 
perform a standard Maxwell construction on the y, vs. 
nE- curves at each gate potential. The resulting phase 
diagram is shown in Fig for back [Fig [^a)] and top 
[Fig [^b)] gate. The dashed orange lines limit the sta¬ 
bility region (in between the system is phase separated), 
corresponding to the densities rq and 712 . We also plot 
the gate voltage as a function of (blue solid lines) for 
different samples having different densities at zero gating. 

As can be seen in Fig. [^the ’’miscibility gap” remains 
open at very negative gate voltage, especially in the back- 
voltage configuration. This does not correspond to the 
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experimental situations [? ]. This discrepancy arises be¬ 
cause, while we took into account the long-range elec¬ 
trostatic cost of the charge distribution, we did not in¬ 
clude in our model any short range repulsion between the 
countercharges (oxygen vacancies or positive charges left 
by the polarity catastrophe reconstruction). While some 
mild fluctuations of electrons and compensating coun¬ 
tercharges are accepted, when too many electrons would 
segregate carrying along too large densities of counter¬ 
charges, the system becomes very rigid and the charge 
segregation stops. 

To overcome this problem, we added a short-range re¬ 
pulsion to the total energy. Of course a precise estimate 
of the energy involved in these repulsive interactions is 
out of reach for the present calculations and we only im¬ 
plement a phenomenological energy barrier that induces 
a rapid grow in /i(n) when electron densities above those 
experimentally found (with a maximum of 0.5el/u.c., in¬ 
cluding the trapped charges) is reached. Due to this 
additional short-range mechanism, the PS region deter¬ 
mined by the Maxwell construction on the dashed red 
lines of Figs. (b-d) is reduced (black solid lines and 
squares in Fig. and, quite remarkably, it ends with 
a critical point at some critical negative gate value, Vg. 
The resulting phase diagram for back-gating is reported 
in Fig.|^(a). The initial (as-grown) density no is deter¬ 
mined by several specific parameters (like, e.g., the num¬ 
ber of LXO planes). Starting from a given no, the total 
density per u.c., n = no -I- Sn, is then changed by the 
gating following the thin blue trajectories in the phase 
diagram. Inspection of these trajectories reveals a very 
large range of intrinsic densities no leading to PS. Upon 
increasing the negative gating the overall average density 
n decreases and the fraction of the system with lower den¬ 
sity rii increases, in agreement with transport measure¬ 
ments [211 nmn], until the system exits the PS region 
at some negative voltage. A similar behavior is found for 
the top-gating case whose phase diagram is reported in 
Fig-ib). 

— Discussion — From Figs.(a) and (b) one can see 
that there exists a broad range of the intrinsic density 
no yielding a negative compressibility, prevented by PS. 
This PS yields an inhomogeneous 2DEG with associated 
inhomogeneity of the oxygen vacancies and/or electronic 
reconstructions. We carried out a detailed analysis (see 
also the Appendix) to identify the specific mechanisms 
determining this instability finding that it arises from 
two main features peculiar to these oxide interfaces. 

First of all, the electrons at the interface are con¬ 
fined in the STO side where a large dielectric constant 
{cr > Coo ^ 100) efficiently screens the electric fields. 
This allows for the accommodation of a large number of 
electrons (« 10^^ cm“^) on several confined levels. The 
large DOS coming from the contributions of the individ¬ 
ual sub-bands greatly enhances the electron compressibil¬ 
ity and facilitates the instability. This effect is stronger 




FIG. 3. (a) Gate potential versus mobile electron density 

phase diagram for the back gated LXO/STO interface in the 
absence (orange dashed lines) and in the presence (black solid 
lines) of short-range rigidity of the countercharges, rii and n 2 
are determined by the Maxwell construction. The blue thin 
solid lines correspond to no = 0.04, 0.05, 0.06, 0.07el/u.c. at 
Vg = 0. (b) Same as (a) for the top-gating case, with no = 
0.04, 0.05, 0.06el/u.c. 


when dxz.yz sub-bands start to be hlled because they have 
a rather large DOS. The relevance of the dxz,dyz levels 
has already been asserted in Hartree-Fock calculations 
m and seems to be experimentally supported [12] • The 
filling of these levels typically starts when the system en¬ 
ters the PS dome upon increasing Vg, and rapidly leads 
to increasingly more abundant regions with locally higher 
electron density. We speculate that this corresponds to 
the observed increase of high-mobility carriers and onset 
of superconductivity mm- 

Secondly, a larger no must correspond to larger density 
of countercharges on the LXO side, which attracts the 
interface electrons and deepens the confining potential 
well, causing a downward shift of the quantized levels. 

Interestingly, in a rather large range of uq values, 
the system exits the PS dome in the vicinity of the 
critical point located here at ng « 0.0475 el/u.c. and 
Vg « —70 V for back-gating (see blue line trajectories 
in Figj^. This suggests that decreasing Vg the electrons 
in the LXO/STO interface eventually display some sig- 
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natures of critical behavior where superconductivity will 
likely be affected by the strong quantum density fluctu¬ 
ations [43]. 

In conclusion, we identified a mechanism of electron- 
driven PS. While the details of the phase-separated re¬ 
gion also depend on the short-range rigidity of the sys¬ 
tem, the existence and robustness of the PS is on a firm 
ground and can be responsible for the strong inhomo¬ 
geneity observed at LXO/STO interfaces. This mecha¬ 
nism can also cooperate with other intrinsic [38] or extrin¬ 
sic (e.g., defects or domain walls UHlIlll) mechanisms. 

We thank V. Brosco, C. Castellani, R. Raimondi, and 
G. Seibold for stimulating discussions. We acknowledge 
financial support from the Sapienza University Project 
n. C26A115HTN, the CNRS PICS program “S2S”, 
the ANR JCJC “Nano-S02DEG”, the SESAME, OXY- 
MORE and CNano programs from Region IdF. 


APPENDIX 


Numerical calculation of the self-consistent 
electronic well 


A well-tested approach to an eigenvalue problem in 
electronic quantum wells at the interface of two semicon¬ 
ductors is reported in [45] . It consist in the simultaneous 
solution of the coupled Poisson and Schrodinger (in the 
effective mass approximation) equations 


2mz dz^ 

e 


-I- e(l){z) + £i 


dz 


er(A) —(/)(z) 


Ci(z) = 0, I = 1,2, 3,... (4) 
= n{z) (5) 


which gives us both the eigenfunctions and the confining 
potential well. The close analogies between semiconduc¬ 
tors and oxides heterointerfaces allows us to apply this 
method to the LXO/STO potential as well. The bound¬ 
ary conditions for the set can be understood con¬ 

sidering the schematic picture of the interface shown in 
fig- Q- Within this approach, the electric field and the 
displacement field are by symmetry on;y directed along 
the z axis. Therefore henceforth we will indicate with E 
and D only the z components of the fields. 

The energy gap between the LXO and STO conduction 
bands is about 3eU, so that in the low-filling regime it is 
safely to approximate the potential step at the LXO/STO 
interface as an infinite energy barrier. In addition the 
eigenfunctions have to be normalizable, therefore we are 
led to set of boundary conditions 


0(0) =0 0(oo)=0. (6) 



FIG. 4. Schematic view of the back-gating configuration in 
LAO/STO. 


The second condition can be implemented numerically 
exploiting the fact we are dealing with a confining poten¬ 
tial: in a region far enough from the interface the eigen¬ 
functions are exponentially decreasing. If Lg resides in 
that region, the set ®> can be replaced by 


0(0) = 0 0(is)=0. (7) 

The boundary conditions for the Poisson equation de¬ 
pend on both the (positive) countercharges left in the 
LXO valence band and on the gating voltage. Experi¬ 
mental evidences show that the width of the 2DEG is 
only few nanometers, while the trapped charges extend 
over a few tens of nanometers mm- If Lp is the thick¬ 
ness of the region in which all the charge is confinedwe 
can solve the Poisson equation analytically in the interval 
(Lp, Ltot) and find appropriate boundary conditions by 
the reasoning reported below. We take Lp = lOOnm, 
that is large enough to study also long-tail distributions 
of the trapped charges. 

Let us consider the interval [0“,-|-oo): in this region 
the system is neutral and we can reasonably impose the 
condition 


E{—oo) = E{+oo) = 0 (8) 

on the electric field, for any gating potential Vq. In 
[0+, -boo) the neutrality does no longer hold, because the 
positive charges (left at the interface of LXO by the po¬ 
lar catastrophe or homogeneously distributed in the LXO 
layer if they are given by the oxygen vacancies) are now 
excluded. Thus the electric field verifies the condition 
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/ dz-^ierE) = +'^^eri0+)E{0+) = -'^ (9) 

Jo+ ^0 Co 

where ng are the polar catastrophe charges and we have 
used Eq. ([^ to cancel the term er(+oo)£'(+oo). Since 
all the electronic charges (mobile n"* = n™ + Sn, and 
trapped rig) are between 0+ and Lp, the electric field 
in z S {Lp, Ltot) is uniform even for the non linear 
experimental form of the relative dielectric constant e^. = 
{A + B\E\)~^ + Coo B and Coo are experimentally 
measured constants mi The boundary value for the 
potential thus reads 

0(Lp) = -^Lp, (10) 

LtOT 

where Vg is the external gating. 

Numerical solution 


that determine univocally the distribution of the 
trapped charges, but a reasonable choice is a de¬ 
creasing exponential distribution of the form 

«o(^) = (y) (11) 

Its shape is controlled by two parameters: the de¬ 
cay length A and the maximum value Ug/A. 

• The Poisson equation is solved in the interval 
[0+,Lp]. This gives a new potential (f){z) = 
4>hew{^) Slid a new electric field E{z) = E^^^{z). 

• The update at the next step is performed defining 
the quantities 

(/)u+l(z) = S(j){^ldiz) + (1 - s)cl)^ewiz) 

E„+i(z) = sE^Uz) + (1 - s)E:,^{z) 


Because of the nonlinear behavior of Q-([^, an ana¬ 
lytic approach is not viable and we adopted a numerical 
iterative procedure to performs the calculation through 
the following steps: 

• The starting point is the potential at the nth it¬ 
eration, and his derivative, the 

electric field En{z) = E";^(z), z e [0, Ls\. We solve 
Eq. Q using a Finite Element Method {FEM) al¬ 
gorithm [46) . finding the energies and the eigen¬ 
functions Ci which verify the conditions Q. We 
recall that Ls has to be taken far enough from 
the interface to ensure that the eigenfunctions are 
in a region in which they have an evanescent be¬ 
havior. This can be done choosing different values 
for Ls and checking them a posteriori. The enve¬ 
lope functions are normalized imposing the condi¬ 
tion \C{z)\^dz = l. 

• From the knowledge of the eigenvalues, the 

Fermi level can be calculated inverting numer¬ 
ically (e.g, using a bisection method) the re¬ 
lation n™ = Jll^g{e)de, and from n"^{z) = 

Ulo d£gi{s) we have the density profile 
along z of the 2DEG. Here only the mobile elec¬ 
trons contribute to define the Fermi energy since 
the trapped charges reside in impurity levels, deep 
or localized enough to be safely considered as a sta¬ 
ble quenched charge distribution. 

• The profile of the potential well is found once one 
knows both the distribution of the 2DEG and the 
trapped (negative) charges. While the first is cal¬ 
culated self-consistently, the second is fixed at the 
zeioth step. There are no experimental evidences 


where the parameter s € [0,1] is used to control 
the difference between the old potential and the 
new one; this reduces the risk that the iterative 
procedure escapes from the self-consistent solution 
and does not converge. The calculation stops as 
soon as the condition 

Lp 

\4>n+l- (ln?dz < 5 ( 12 ) 

is fulfilled, for a suitably chosen accuracy <5. 


Through the previous five steps we obtain the poten¬ 
tial well V{z) = —e(j){z) (with the corresponding Fermi 
energy cp), the electric field E{z) and the density nE‘{z) 
of the 2DEG for a given value of ng. 

In principle the numerical solution may depend on the 
choice of the error 5 in (12), the initial potential 4>o{z) 
and the discretization length (according to the standard 
FEM |46], we discretized in Neiem intervals the region in 
which we solve the Schrodinger equation). We performed 
different tests to ensure the reliability of the numerical 
solution. In order to exclude the dependence on the ini¬ 
tial potential, we varied the initial condition (jjQ in a rea¬ 
sonable class of functions and kept ng and Vg fixed. We 
find that the numerical solution is stable with respect to 
the initial data. IVejem and 6 has been fixed to reach 
the uniform convergence. The solution is almost inde¬ 
pendent of 6 while the choice of the first parameter is 
critical to obtain the uniform convergence. In our work 
we took Nelem SO that sup„j, \eF{Nelem + 1) “ eF{Nelem)\ 
is of order « 10“®, a suitable request to reduce the error 
under the typical variations (~ 10“^eE) of the chemical 
potential /r(n). 
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Electrostatic energy 


While the previous section solved the quantum prob¬ 
lem of electrons in a self-consistent potential, here we pro¬ 
vide details of the calculation of the other electrostatic 
contributions to the total energy of the system arising 
form the fixed charges (those on the gating electrodes 
and the charges trapped in the impurity states inside the 
STO) and from the mobile charges at the interface. 

We consider two possible gating configurations. For 
the back-gating, the near region is the interval z G [0, Lp], 
while the far region is the z G (Lp, Ltot] interval. On the 
other hand, for the top-gating, while the so-called near 
region is the same, the “far” region corresponds to the 
LXO side of the heterostructure and corresponds to the 
z G [—d, 0) interval. 

In the near region we solve the Poisson equation 


^Diz) = n"^{z)+nl(z) (13) 

to find the electric displacement field D{z) = D^{z) + 
Df{z) due to the mobile and fixed charges. lin this re¬ 
gion the only fixed charges are those trapped in the STO 
side, while the charges on the electrodes and the coun¬ 
tercharges on the LXO side only enter as boundary con¬ 
ditions. The electric field is 


E = 


+ 


(14) 


- sgn(£>) [1 -t eppA - sgn{D)BD] 

‘ZBcqCqq 

-I- Coo A - sgii{D)BDf + sgii{D)4eooABD 

271eo6oo 


The total electric displacement field is obtained by nu¬ 
merically integrating Eq.(13), while the field due to the 
fixed charges is analytically found from 

d 


dz 


D^{z) = no(z). 


(15) 


The result is 


Df (z) = - 


—en\ 


n / 2 \ ^^0 gStI . 

+ (M) 


Once D and D-^ are known, the electric displacement 
due to the mobile charges is obtained from D™ = D—D^. 
The electric fields of the fixed and mobile charges are 


£)/,™ A + B\E\ 

Co 1 + Coo (^ + -B|£'|) 


(17) 


Notice that the non-linearity of Cr entails the dependence 
of D-l’™ and on the total electric field E given by 


Eq.(|^ 

In the far region we have different expressions for the 
top- and back-gating configurations. Eor the latter the 


far region coincides with the bulk of the STO substrate, 
where 


£i'"(z) = 
D^{z) = — 


e»o 

2a2 

e»o 

2a2 


eSn eno 
eno edn 


(18) 

(19) 


are constant. The corresponding electric fields are then 
obtained by solving Eqs.(14) and Q- 

The far region in the top-gating configuration is in¬ 
stead given by the whole LXO side, where the dielectric 
constant eLxo is really constant 


E™(z) 

Efiz) 


edn 

‘2-eoeLXOOp' 

edn 

‘^eoeLXoa^ 


e»o 

2eoeLXoa^ 

e^o 

2eoeLXoa^ 


eno 

2eoeLXoa^ 

enp 

2eoeLXoa^ 


( 20 ) 
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